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ABSTRACT 

All atmosphere-less planetary bodies are covered with a dust layer, the so-called regolith, which 
determines the optical, mechanical and thermal properties of their surface. These properties depend 
on the regolith material, the size distribution of the particles it consists of, and the porosity to which 
these particles are packed. We performed experiments in parabolic flights to determine the gravity 
dependency of the packing density of regolith for solid-particle sizes of 60 /rm and 1 mm as well as for 
100-250 /rm-sized agglomerates of 1.5 /rm-sized solid grains. We utilized g-levels between 0.7 m s“^ 
and 18 m s“^ and completed our measurements with experiments under normal gravity conditions. 
Based on previous experimental and theoretical literature and supported by our new experiments, we 
developed an analytical model to calculate the regolith stratification of celestial rocky and icy bodies 
and estimated the mechanical yields of the regolith under the weight of an astronaut and a spacecraft 
resting on these objects. 


1. INTRODUCTION 


Gary et al. (1972) state that “Regolith is a mantle of 
loose incoherent rocky material of various origins that 
forms the surface of planetary bodies”. Its size distri¬ 
bution determines the optical, mechanical and thermal 
properties of these surfaces. On bodies without atmo¬ 
sphere, regolith is formed by high-velocity impacts of in¬ 
terplanetary particles of different sizes. These impacts 
produce ejecta material with a size and velocity distri¬ 
bution determined by the impactor and target proper¬ 
ties. Ejecta with velocities lower than the escape ve¬ 
locity of the target body are reaccreted and, thus, form 


Hartmann 

(1985 

1 found ejecta velocities depending on 
le impactor. The experiments of (e.g. 

the velocity ot tl 

Fujiwara and Tsukamoto 

1980 

Nakamura et al. 

1994) 


and studies of crater structures on the Moon (V ickery 


1986 1987) showed that smaller particles are ejected at 


higheiUleloc ities than the larger ejecta. Finally |Bottke| 


et al. (2002) calculated in their chapter 111 the trajec¬ 
tories and re-accretion of ejected material on small ce¬ 
lestial objects depending on its velocity. Therefore, the 
regolith size distribution should depend on the escape 
speed and, hence, on t he size of the target body. |Gund-| 
lach and Blum (2013) developed a method to correlate 
the size of the regolith particles with the thermal con¬ 
ductivity of the regolith. Thus, from measurements of 
the thermal inertia of small objects in the Solar System, 
their regolith-particle size can be estimated. The regolith 
filling factor, i.e. the packing density of the regolith par¬ 
ticles (or 1 - porosity), enters the i r calc ulation as a free 


parameter. Gundlach and Blum (2013) found that ob¬ 


jects with a size smaller than 100 km are covered with 
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regolith particles of the size of ~1-10 mm, whereas larger 
o bjects carry regolith with particle sizes of 10-100 /rm. 


Skorov and Blum (2012) model comet nuclei as con- 
sisting ot m acroscopic dust and ice agglomerates. In a 
recent work, |Blum et al. (2014) showed that this model 


can explain the observed continuous activity of comets 
in the inner Solar System. Inside the comet nucleus, the 
dust and ice agglomerates are only subjected to the weak 
gravitational force of the comet, which leads to a stratih- 
cation of the packing density within a comet nucleus. On 
celestial objects the regolith is compacted by their grav¬ 
itational force, which leads to a stratification where the 
filling factor of the regolith increases with layer depth. 

However, the texture (e.g. particle size distribution 
and packing density) and vertical stratification of re¬ 
golith is basic ally unknown for all celestial ob jects except 
for the Moon (Heiken et al.|[l9^). Recently, [Kiuchi and 


Nakamura (|20r3) related the particle size to the porosity 
for a variety ot celestial objects. However, they used a 
simplified model in their analysis that implies that the 
packing density depends only on the ratio of the gravi¬ 
tational force of the contact force of a regolith particle. 
A result of this simplihcation is that the packing density 
of the regolith does not change with depth. If this was 
true, it would be impossible to co mpress loosely packed 
regolit h at all. Thus, the results of |Kiuchi and Nakamura] 
(2014) are valid for the uppermost layer of the regolith 
where hydrostatic compression is small. 

We will show in this article that stratification of re¬ 
golith due to hydrostatic compression is important and 
can be calculated with the relation between pressure and 
filling factor of the regolith found in our article. With 
this information, also the yield of objects on the surface 
of the regolith-covered celestial objects was derived. 

For /rm sized grains, this relation is experimentally ob¬ 
tained (see Section]^ by compressing the regolith with a 
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piston in a cylinder as shown by Guttler et al. (2009). For 
the 10 fim grains, we measured the low-pressure part of 
the compression curve by adding thin layers on a regolith 
to compress it by its own weight and the high pressure 
part by compressing the regolith by a piston. For grains 
larger than 50 p,m, the regolith is compacted to RCP 
(random close packing, its highest possible filling factor) 
by its own gravitational force for thicknesses above a few 
particle layers. We therefore did most of these measure¬ 
ments in the low-gravity environment aboard the Zero-G- 
plane of the European Space Agency. However, in spite 
of that, it was only possible to meas ure the upper 20 % 
of the compression curve (see Section 2.31. It was there¬ 
fore necessary to develcm a model for the full compres¬ 
sion curve (see Section ^ in which only one parameter 
must be experimentally aetermined. With this analyti¬ 
cal model, we could additionally reproduce the re lation 
between fill ing factor and grain radius found by |Yang| 
et al. (20001 at constant and can show the consistency 
of our model (see Sect. [^. 

To simulate the stratmcation of the packing density 
within a comet nucleus, we additionally used in our ex¬ 
periments 0.1-0.25 mm-sized dust agglomerates consist¬ 
ing of 1.5 /rm-sized mono-disperse and spherical Si 02 
grains. We performed our measurements with mono- 
disperse spherical monomer particles to ease the theoret¬ 
ical approach to our measurements, whereas the particles 
on celestial objects will most likely be irregular and poly- 
disperse. H owever, our previous measurements in |Blum| 
et al.| (|2006D showed that the results wil]_not change dra¬ 
matically tor these particles (see Sect. 

In Section we describe our experimental approach 
and the experimental findings. Section uses the 
analytical-approx imation form of Gutt l er et al. ( 2009|) 
and the resu lts of |Dominik and Tikens 1 1995) and Krijt 
et al. (2014) to rnodel the observed regolith stratihca- 
tion. In Section]^ we derive filling-factor profiles of the 
regolith for selected small rocky and icy bodies and es¬ 
timate the mechanical yield of the regolith under the 
weight of an astronaut and a spacecraft, resting on these 
objects. Finally, Section summarizes our results. 

2. EXPERIMENTAL APPROACH 


In this Section, we describe the experimental meth¬ 
ods developed and applied for the determination of the 
pressure-dependent filling factor of regolith and the re¬ 
sults we received in parabolic-flight and laboratory ex¬ 
periments. 


2 . 1 . Particle Samples 

As regolith analogs, we used amorphous Si 02 spheres 
with diameter of 10 /rm, glass spheres with 60 fim diam¬ 
eter from unknown glass type, spherical soda lime par¬ 
ticles with 1 mm diameter, and agglomerates of 1.5 /rm- 
sized amorphous Si 02 spheres with 100-250 ^m diameter. 
Fig. shows the size distributions and microscopic im¬ 
ages of these samples. The experiments were performed 
at 10® Pa ambient pressure. Laboratory, microgravity 
and hypergravity experiments were performed with all 
samples with the exception of the 10 ^m-sized particles. 
These particles fluidized under microgravity conditions, 
due to the ambient gas pressure, and could, thus, not 
deliver reliable results. For these particles, we only ana¬ 
lyzed the experiments conducted in the laboratory. The 
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Fig. 1.— Measured number-fraction histograms of the diameters 
of the regolith-analog particles (left column) and the corresponding 
microscopy images of our particle species. The inset in the bottom- 
right picture shows the constituent 1.5 /rm-sized Si02 grains of the 
dust agglomerates. 

humidity of the air inside the experiment could increase 
the contact forces of the particles and therefore decrease 
their filling factor at a given pressure. To find this in¬ 
fluence we did ground experiments for the 10 /im-sized 
particles at 10 Pa gas pressure and at ambient (10® Pa) 
air pressure (asterisks and triangles in Fig. and found 
no deviation within our error bars. 

2 . 2 . Experimental Setup 
2.2.1. Parabolic-Flight Experiments 

Parabolic-flight experiments under reduced and hyper¬ 
gravity conditions were performed onboard the ZERO- 
G Airbus A300 aircraft. Typically, several experimental 
setups are mounted in the passenger area of the aircraft 
and are operated by onboard experimentalists during the 
flights. In one such flight, the aircraft performs up to 31 
parabolic and catenary-curve flight maneuvers. 

During a parabolic flight maneuver, inertia forces com¬ 
pletely cancel Earth’s gravity and the aircraft and all ex¬ 
periments onboard are completely weightless. During a 
catenary-curve flight maneuver, inertia forces partly can¬ 
cel Earth’s gravity and the aircraft and all experiments 
are partly weightless so that the environment of celestial 
bodies with surface accelerations smaller than on Earth 
can be simulated. 

In the flights described here, the flight maneuvers con¬ 
sisted of 13 catenary curves with Martian (0.38 5 ), 12 
with Lunar (0.17 g) and 6 parabolas at zero gravity, 
respectively. Here g = 9.81 m s“^ is the surface ac¬ 
celeration on Earth. As the aircraft has to enter and 
exit the parabolas and catenary curves, a ’’pull-up” 
and a ’’pull-out” maneuver are required before and af¬ 
ter each constant-acceleration curve. These maneuvers 
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cause hyper-gravity accelerations from 1.3 g to 2 g. We 
also utilized these hyper-g phases for our experiments. 

In the parabolic-flight experiments considered here, 
four different masses of otherwise identical regolith sam¬ 
ples were filled in glass cylinders with an inner diameter 
of 2.5 cm, with cylinder no. 1 having the highest mass 
and cylinder no. 4 with the lowest mass, respectively. 

Even for the largest (1 mm) particles, only 1% of their 
contact points are with cylinder walls so that the finite 
cylinder diameter will infl uence the measurements only 
through the Janssen effect ( Jansen[1895| 
tihed m Eq. It 


Sperll2005l de- 


should also be 


scribed below and quantihea 
noted that the lowermost layer of particles naturally pos¬ 
sesses a different filling factor, because of the flat surface 
of the cylinder base. This boundary effect is relevant for 
the mm-size particle measurements and was corrected for 
in the data analysis. 

The four glass cylinders were placed in a polycarbonate 
experimental box and mounted on a linear stage together 
with a video camera (see Fig. 1^. The video camera was 
used to measure the filling height of each of the samples. 
Caused by the mounting of the cylinders in the polycar¬ 
bonate box, their bottom was not visible in the field of 
view of the camera (see Fig. [^. Therefore, we calibrated 
the hlling-height by comparing the measurements during 
the 1 g phases with height measurements on the ground 
without the box. 

During the parabolic flights, we performed two types 
of experiments. 


• Type-I experiments: The regolith is deagglom- 
erated and homogeneously distributed inside the 
glass cylinder by a single shake using a succession 
of an upward 20 g acceleration, a downward -20 g 
acceleration, and then again an upward 20 g accel¬ 
eration so that in the end the sample box comes 
to rest at its starting position. This initial shaking 
with an amplitude of 20 cm was performed using 
the electromagnetic linear stage. Thereafter, the 
regolith particles sedimented in the ambient (log-g 
or hyper-g) acceleration downwards and formed a 
layer whose filling height could be measured on the 
video-camera images (see Fig. [^. During this sed¬ 
imentation phase, the residuaFacceleration vector 
in the aircraft is not perfectly directed in vertical 
direction. Therefore, surfaces of the regolith in the 
4 cylinders are slightly tilted (in all 4 cylinders in 
the approximate same direction and with the same 
angle). We only used measurements in which the 
surface was tilted towards the camera so that the 
full perimeter of the sample surface was visible. We 
measured the fill height at four points at the cylin¬ 
der rim, that have a distance of 90 degrees and 
calculated the mean fill height by averaging over 
these four measurements. A random deviation of 
the surface from a perfect plane is a reason for the 
large error bars in our figures andbecause an 
error in the fill height measurements of ±0.3 mm 
corresponds to an error in $ of ±0.1. The error 
bars of this measurements are found statistically 
by calculating the mean and standard deviation of 
different measurements. 

The ambient acceleration level is produced by the 
airplane that performs catenary-curves with Mar¬ 


tian, Lunar and hyper-gravity accelerations. As the 
acceleration level provided by the aircraft varied 
strongly with time, particularly during the hyper¬ 
gravity phases, and as the experiment duration 
(i.e., shaking plus sedimentation) was always short 
compared to these variations, we were able to mea¬ 
sure the fill heights for about 50 uniformly dis¬ 
tributed g-levels between ~ 0.3 g and ~ 2 g during 
the parabolic-flight campaign. For ease of analysis, 
we equipped our experiment with an acceleration 
sensor whose read-out value (in units of g) was dis¬ 
played in the field of view of the video camera (see 
the LED display in Fig. [^. 

• Type-2 experiments: These experiments were per¬ 
formed at micro-gravity conditions generated by 
parabola maneuvers of the airplane. The regolith is 
deagglomerated and homogeneously distributed in¬ 
side the glass cylinder as described in the type-1 ex¬ 
periments. In a subsequent second step, the linear 
stage is very slowly (at 5 cm s“^) moved upwards 
to bring all regolith particles back to the base of the 
cylinder. In a third step, the linear stage acceler¬ 
ated the experimental box upwards to compress the 
granular medium at pre-defined constant accelera¬ 
tion levels of 0.05, 0.1, 0.15 and 0.2 g, respectively, 
over a length of 71 cm, causing acceleration times 
of 3.7,3.1 and 2.7 s, respectively. Naturally, these 
accelerations were superposed with the residual ac¬ 
celerations of the aircraft, which were in the range 
of ± 0.2 g. We were able to perform up to three 
type-2 experiments per parabola. Because the tem¬ 
poral changes of the residual acceleration were long 
compared to the experimental duration, we again 
obtained a total of around 40 equally distributed 
g-levels between ~ 0.06 g and ~ 0.3 g. 


2.2.2. Ground Experiments 

For the experiments exclusively performed in the lab¬ 
oratory under normal-gravity conditions, we successively 
filled a small mass dm of regolith simulant particles in a 
glass tube with 25 mm diameter and measured the result¬ 
ing gain in volume, dV. The ratio dm/dV is the density 
of the lowermost level in the glass tube. Repeating this 
procedure by carefully stratifying additional layers with 
a spatula resulted in a filling-factor proHle along the tube 
height (type 1 ground experiments). 

As it was not possible to perform reduced-gravity 
experiments for the 10 gm sized particles during the 
parabolic flights, we additionally conducted type 2 
ground experiments. The difference to the regular (type 
1 ) ground experiments was that we additionally shook 
the regolith before measuring the fill height to avoid any 
hysteresis effect. To obtain high pressure data of the 
10 gm sample, we al so performed compr ession measure¬ 
ments as described in|Cuttler et al.|(|2009|) (type 3 ground 
experiments). We performed these experiments in vac¬ 
uum at 10 Pa gas pressure (triangles in Fig. and at 
ambient air pressure (lO^Pa; asterisks in Fig. ^ and 
found that there is no difference within our error oars. 

An advantage of using 10 gm sized particles over larger 
grains is that the hydrostatic compression by their own 
weight is much smaller than the pressure required to com- 
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Side View Front View 



Camera View 



Fig. 2. — Top: Schematic drawing of the experimental setup. 
Bottom: Picture of the experimental box with the four sample 
cylinders containing regolith from 60 pm-sized spheres, captured 
with the video camera. The LED display on the upper left of the 
image shows the acceleration level inside the airplane. 


press them to RCP. Therefore, higher hydrostatic pres¬ 
sures were simulated by compressing the regolith with 
a piston. Thus, we could achieve measurements over a 
wide range of the compression curve for the 10 /rm sized 
particles. For the 60 fim and 1 mm particles, only seg¬ 
ments of this curve could be obtained (see Fig. |3. 

The pressure acting within the regolith bottornlayer in 
the cylinder scales with the fi ll height x of the cylinder s 
following Janssen’s Equation ( Jansen|fl895 Sperl|[2005 ) 


y n 

P - [ 


exp(-4iV^)], 


( 1 ) 


where K = 0.2 (Jansen 1895| ), TO, and r are the ratio 
between the horizontal and vertical components of the 
stress tensor in the regolith, the mass inside the cylinder 
at a given fill he ight x, a nd th e radius of the cylinder, 
respectively. As Jansen (18951 used containers with a 
square cross section but we employed circular cylinders 
instead, we substituted the side length by 2 r. 


2.3. Experimental Results 
2.3.1. Parabolic-Flight Experiments 


As shown in Sect. 2 . 2 . 1 [ our measurements in 
the parabolic-flight experiments directly yielded the g- 
dependence of the average filling factor of the regolith in 
each of the four glass tubes. As these four glass tubes 
were always filled with the same sample material but to 
differ ent heights, we used the method depicted in Sect. 
|2.2.2| to derive the height-dependent filling factor. Thus, 
cylinder 4 with the smallest fill height corresponds to the 
uppermost layer of the other three cylinders and pos¬ 
sesses an average filling factor of 


$4 = 


TO4 

Lpb’ 


( 2 ) 


where m 4 , V 4 , and pf, are the mass and volume of the re¬ 
golith in cylinder 4, and the density of the regolith bulk 
material, respectively. The filling factor of the subse¬ 
quent layers f = 3, 2,1 can then be calculated by 


TOi - TOi+i 

(V-V+i)pb’ 


( 3 ) 


with TOi and Vi being the mass and volume of cylinder i, 
respectiv^. 

Figs. shows the resulting filling factors (i = 
4,3, 2,1 from top to bottom) of the regolith layers con¬ 
sisting of 60 /im-sized particles, 1 mm-sized particles, and 
100-250 pm-sized dust agglomerates, respectively, as a 
function of the ambient acceleration level. The error bars 
denote the standard derivation among the 10-12 (60 pm- 
sized particles), 16-20 (1 mm-sized particles), and 17-21 
measurements (100-250 pm-sized dust agglomerates), re¬ 
spectively. 

The graphs for the 60 pm particles generally show 
an increase of the filling factor with increasing depth 
and increasing gravity level, until RCP (<i)RCP=0.64) is 
reached. The filling factor for the mm-sized particles is 
close to RCP for the whole gravity and depth ranges. As 
the individual measurement errors are very large, it is 
difficult to recognize a systematic behavior for these par¬ 
ticles. The interpretation for this data is easier in Fig. 
1 ^ because there the data of layer depth and gravity level 
are combined to give a better statistical basis. The dust 
agglomerates are, however, only compressed to a max¬ 
imum filling factor of $ = 0.37 and never get close to 
RCP. 

In all our measurements, the regolith particles adhered 
in mono-layers to the walls and to the top of the cylin¬ 
ders. However, this considerably affects the resulting 
filling factors only for the mm-sized particles. We could 
partly correct for this effect for the 1 mm-sized spheres by 
counting the (spatially resolved) particles on the walls, 
but not all surfaces were visible in the camera field of 
view. For the invisible surfaces, we assumed the same 
surface coverage of adhered particles as on the visible 
surfaces. However, we have to bear in mind that this 
procedure might lead to an overcorrection of the filling 
factor, because the invisible surfaces were in the upper 
part of the cylinders and, thus, not necessarily came in 
contact with the same number of particles. 

The effect of adhesion is stronger at low p-levels and 
leads to the artificial decrease (because of the overcorrec¬ 
tion described above) of the filling factor for accelerations 
between 0.1 g and 0.3 g (see Fig. 1^, which is visible in 
the data of the 1 mm-sized partiaes. These particles 
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also show a boundary effect in cylinder 4, because of the 
rather small layer thickness of 4 mm and a naturally re¬ 
duced filling factor in the first particle layer from the 
cylinder bottom. 

It was not possible to use 10 /im sized particles in the 
parabolic-flight experiments. The inside of the cylinders 
were at ambient pressure and an air cushion formed in 
between the particles. The time available for sedimenta¬ 
tion during the experiment was too short for the air to 
diffuse out of the regolith layers. This effect also led to 
a systematic reduction of the Ailing factor for the 60 /rm 
particle measurements at low pressures. 

However, the model described in Section will show 
that the Ailing factor only depends on the hydrodynamic 
pressure inside the regolith. Therefore, in Fig. we 
plotted all reduced- and hyper-gravity data points as a 
function of the respective hydrostatic pressure inside the 
regolith, using Fq. Q, togeth er wit h the data from the 
ground experiments ^ee Sect. 2.3.2). 


2.3.2. Ground Experiment 

Microscopically, restructuring and compaction inside 
a regolith layer takes place, because external forces act 
against rolling fric tion on inter-particle contact points 
(Yang et al. 2000). These forces arise from the ambi¬ 
ent pressure inside the regolith and their origin (gravity- 
driven hydrostatic or external by a piston) is irrelevant. 
Therefore, the compression of regolith can also be studied 
with ground-based experiments for pressures exceeding 
the hydrostatic value. The results of our type 1 ground 
experiments, as described in Sect. 12.2.2 are shown as 
pluses in Fig. and were derived by boxcar-averaging 
of four individual filling-factor measurements as a func¬ 
tion of hydrostatic pressure. We chose to use a boxcar¬ 
averaging method instead of binning the data points 
(which was done for the parabolic-flight experiments), 
because the experiments showed a hysteresis effect, which 
would have been veiled by the binning method. Be¬ 
tween two subsequent measurements in the laboratory, 
we a dded a defined mass of particles to the sample (see 
Sect. |2.2.2|). 

In the type 2 ground experiments measurements, the 
regolith was shaken up after a new upper layer had been 
deposited so that all regolith particles could sediment 
simultaneously. With this, we got rid of the above- 
described hysteresis effect and could bin the data for 
better statistics. These data (only available for the 10 

m particles) are shown in Fig. |^as triangles (measure¬ 
ments in vacuum) and asterisks (measurements at 10® Pa 
air) with error bars. In the type 3 ground experiments, 
we produced high pressures by compressing the re golith 
with a piston as described in Guttler et al.| (|2009|) . The 
results of these experiments (only tor 10 /r m-particles) 
are shown as pluses with error bars in Fig. Here, each 
error bar denotes the standard deviation of rour different 
measurements. These measurements start at around 10^ 
Pa pressure. As can be seen in Fig. slightly increasing 
pressure does not appreciably compact the regolith. We 
think that is because the internal friction of our compres¬ 
sion apparatus has to be overcome. At around 10® Pa 
pressure, the regolith then gets compacted. At 10® Pa 
pressure, our measurements terminate, because the limit 
of our experimental apparatus had been reached. 

Fig. I^shows a general increase of the filling factor with 
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Fig. 3.— The filling factor $ for regolith layers consisting of 
60 /.fm-sized spheres (left), 1 mm-sized spheres (middle), and 100- 
250 /im-sized dust agglomerates (right), respectively, as a function 
of the ambient acceleration level (in units oi g = 9.8 m s“^) at 
depths of 3 mm, 7 mm, 10 mm and, 14 mm (from top to bottom; 
left), 2 mm, 5 mm, 9 mm and, 13 mm (middle), and 3 mm, 6 mm, 
7 mm and, 10 mm (right), respectively. Each data point in the 
graph represents the mean of 10 - 12 (left), 16 - 20 (middle), and 
17 - 21 (right) individual measurements, respectively, and the error 
bars denote one standard derivation. 
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increasing pressure, but in some intervals of the type 1 
ground experiments, the filling factor decreases with in¬ 
creasing pressure. This is due to the fact that a granular 
medium does not always restructure under the weight of 
the additional layers, which then leads to a reduction in 
filling factor. After the additional layers reach a height 
between 4 mm and 7 mm, corresponding to an additional 
pressure of 20 to 35 Pa, the regolith suddenly restructures 
to a higher filling factor. The effect is exaggerated in the 
data shown in Fig. because the hysteresis, which takes 
place over the whole cylinder height, is assigned to only 
one layer, due to our analysis method. Further evidence 
to this behavior is visible in the data, which do not show 
arbitrarily scattered data points but follow an undulat¬ 
ing line, which is clearly visible for the agglomerate data 
at high pressures (pluses in the rightmost panel of Fig. 
1^. For the other particle types, this behavior is only 
visible when the data points are plotted linearly. 
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Fig. 4.— The filling factor <I> for regolith layers consisting of (from left to right) 10 /rm-sized, 60 /rm-sized, and 1 mm-sized spheres as 
well as for 100-250 /rm-sized dust agglomerates as a function of the hydrostatic pressure, following Eq. This pressure was calculated from 
the layer depth and the ambient acceleration level as measured in the experiments. Triangles represent data points from the uppermost 
layer, diamonds from the second layer, asterisks form the third layer, and crosses from the bottom layer (for the 60 fim and 1 mm-sized 
spheres and for the 100-250 /rm-sized dust agglomerates only). These data points and their error bars were measured in the parabolic-flight 
experiments and are identical to those shown in Fig. 1^ Pluses represent boxcar averages of the filling factor measured in the type 1 ground 
experiments. Mind that for the 10 /rm-sized particle, only ground-based data are available. For these particles, the data measured in 
type 2 ground experiments are shown by triangles (vacuum) an asterisks (air) and for the type 3 ground experiments by pluses with error 
bars. The solid line corresponds to the model curve as given by Eq. It is fitted to our measurements using pm as a fit parameter. The 
upper and lower dashed lines correspond to Eq. |^with Pmin Pmax, respectively. Mind the different pressure scale for the 10 /rm-sized 
particles. 


Another issue of the type 1 ground experiments is that 
after the addition of another particle layer, the sample 
had to be moved to determine the new mass. Any as¬ 
sociated vibrations may have caused restructuring and, 
thus, may have artificially increased the filling factor of 
the sample. The data of the mm-sized particles are again 
affected by the loss of particles, which stick to the walls 
at low and medium pressures, which increases the mea¬ 
sured filling factor. Apart from this, our laboratory data 
fit well to our parabolic-flight results. 


3. MODELING 

3.1. Regolith consisting of solid spherical particles 


Guttler et al. (2009) used the following empirical de- 
scription tor the pressure dependence of the filling factor 
of cohesive-particle layers: 

$(E) = $2 - 


$2 - $1 


exp( ‘°g^^°gP- ) 


1 


(4) 


Here, S is the pressure and $i, $2, Pm, and A are the 
minimum and maximum filling factor at very low and 
very high pressures, the turnover pressure and the log¬ 
arithmic width of the transition from low to high fill¬ 
ing facto r. The trans i tion w idth was experimentally 
found by Guttler et al. (2009) to be A = 0.58 for omni- 
directional compression. I’he packing density $1 (Ran- 
dom Ballistic Deposition, R BD) has been measured by 
Blum and Schrapler (2004) which matches results from 
numerical simulations and corresponds to a coordination 
number of 2 for spherical monodisperse particles, inde¬ 
pendent of their radius. RBD means that particles se¬ 
quentially drop at random positions on a surface and 
stick where they hit. $2 (Random Close Packing, RCP) 
is used for spherical monodisperse particles independent 
of their size and corresponds to a mean coordination 
number of 6 (Yang et al. 2000). At RCP, a regolith 
is jammed and cannot be compressed further. 


Sliding is not possible either, because sliding forces are 
very high and the regolith particles will crunch prior to 
sliding. 

At a filling factor of 4) = 0.15, the average co ordination 
number of t he particles inside the regolith is 2 (Lagemaat 
et al. 2001) so that the particles are free to rnove. Any 
pressure is transformed into a force at the particle con¬ 
tacts, which is proportional to the particle cross section 
and the inverse filling factor. With this consideration, 
the pressure at which the restructuring of the regolith 
starts, can be calculated by 


Pi = 


Froll^l 


(5) 


with 


Froll — -^0 



( 6 ) 


(Krijt et al. 2014). Here, is the particle radius. 

Equation 5] relates the pressure inside the regolith to 
the force exerted on a single particle inside the regolith 
that makes it roll over its contact points. To derive this 
force, we multiplied the ambient pressure inside the re¬ 
golith by an effective cross section per unit monomer 
particle. This effective cross section (at a coordination 
number of 2) is the cross section of a monomer parti¬ 
cle divided by the packing density, i.e. the smaller the 
packing density the higher the force on a single particle 
contact for a given pressure. In other words, pi is the 
maximal pressure a regolith can sustain without restruc¬ 
turing. The material paramete r Fo at the radius rq are 
taken from the measurement of Heim et al. (1999). With 
the assumption that A is independent of particle size, 
the ratio 



( 7 ) 
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is particle size independent and it is possible to calculate 
Pm from Pi for all particle sizes, 


Pm — 


Froll^l 


R. 


( 8 ) 


The measurements by Guttler et al. (20091 show filling 
factors slightly below d) = 0.15 at low pressures, which 
implies that their regolith at low pressures included void 
spaces. Therefore, it is not possible to use pi from Eq. 

in this case. For pi we instead have to use the filling 
factor d> = 0.15 to assure that the coordination number 
of the particles is 2. Th erefore we take the p ressure from 
the compression curve of Guttler et al. (2009) at $ = 0.15 
and associate it with pi. 

We assume that the logarithmic transition width A 
is independent of particle size so that we obtain a con¬ 
stant value of A = 0.58 and a constant R for all sam¬ 
ples. In a descriptive way, this is the case, because the 
increasing resistivity against compression at increasing 
filling factors is due to a reordering of the particles and 
a corresponding increase of the mean coordination num¬ 
ber. The mean coordination number is independent of 
the particles size and is directly connected to the force 
chains in a regolith, which are responsible for its resis- 
tan ce to compression. For a proof please compare figure 
4 of Langemaat et al. (2001), which was done for mono- 
disper s e nan ometer sized particles with Figure 9 of|Yang| 


et al.| (|2000|) which is plotted for different sizes (1 /rm 
to 100(J ^mj depending on porosity. These curves are 
nearly identical. 

To prove the assumption of a constant A and to show 
the consistency of our model, we plotted our data for a 
consta nt ambient pressure S (calculated from Yu et al. 
(19971 as described below) ^ a function of the particle 
radius (solid curve in Fig. [^. This was done by eval¬ 
uating the relation between pm and the particle radius 
through combining Eqs. [5l[6]and[3 


Pm{rp) = Tb$irp ^tt Vp 


(9) 

The parameters of the resulting solid curve in Fig. [^are 
shown in Table [T] 

A potential radius dependency of A was realized by 
multiplying A with 


or 


0.5 


-0.5 


( 10 ) 

( 11 ) 


with Tj = 10 pm being a scalin g p arameter. The dotted 
amsnd dashed curves in Figure show the behavior of 
the filling factor as a function ofparticle radius for these 
cases. Both deviations from constant A result in a non¬ 
monotonic behavior of the curves, in contradiction to 
all available measurements and numerical simulations, 
which are monotonically increasing and do not fit these 
curves (see, e.g. the data of sev eral authors given in 
|Yu et alT]|2003[ [Yang et al.||2000[ an d the data shown 
in Figure ]5|f The measurements by Yu et al. (19971 
were performed with irregular particles with a narrow 
size distribution. Their data are shown as pluses in Fig. 


[^. Their experiments and the corresponding setup are 
so well described that we could use their measure ment 
results to calculate the pressure. I Yu et al.l (|1997|) used 


the standard method: they filled a truncated cone with 
an upper diameter of 57 mm, a lower diameter of 44 mm 
and a height of 70 mm to the rim with their particles and 
measured their weight. We calculated the mean pressure 
inside the container for the average mass they filled in 
for different particle sizes and therefore packing densities. 
We did this by integra ting Equation[T]over the fill hei ght. 


Feng and Yu (19981 and Milewski and Katz (19871 did 
measurements with spherical p articles usi n g the same 
method. The measu rements of Yu et al . (|l99 7) (dia¬ 
monds in Fig^. and Milewski and Katz 1 1987|) (trian- 
gles in Fig. fit our calculations properly, which shows 
that our Eq. a constant A, a power law for Pm('Cp) 
as well as the ansatz with a Fermi function is appropri¬ 
ate. The measurements of Yu et al. (1997) show that 
the power law also holds for irregular particles, and that 
irregular particles behave like spherical particles that are 
somewhat smaller. Quantitatively, this means that the 
corresponding pm of irregular grains is approximately a 
factor of 8 larger than for spherical particles so that the 
curves in Fig. are shifted by a factor 5 to the right to 
be valid for irregular particles. With the above consider¬ 
ations, and having in mind that Pm for regolith particles 
with a wide size distribution has the opposite effect on 
Pm and nearly cancels the change from spher ical to ir¬ 
regular grains (see the discussion in Section 3.3), this 
estimates a maximum error of our curves caused by the 

assumption of spherical particles. _ 

The dotte d curve in Fig. [^is a best fit taken from Yu 


et al. (20031 to the measurements with irregular particles 
from several authors. 



Fig. 5. — The dependency of the filling factor $ on the parti¬ 
cle radius at constant fill height and gravity. Pl uses: data from 
irregular particles with a narrow size distribution jYu et al.|1997l|. 
Diamonds: dat a from spherical mono-disperse particles (|Miiew^i| 
land K atz|1987ll. Triangle s: data from spherical mono-disperse par- 
ticles I p'eng and Yu|1998| |. Solid curve: our model with constant A. 
Dashed-dotted curve: our model with A oc r®'^. Dashed line: our 
model with A cc Dotted line: best fit of measurements wi th 

irregular particles of several authors taken from [Yu et al.| | |2003[ l. 


To derive from the measurements shown in Fig. ffl 
we fit Eq. to the measured data, with pm being the 
only fit parameter. The resulting values for pm as well as 
the other parameters used in Eq. [^are shown in the first 
three lines of Table The best-fitting curve is shown as 
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a solid line in Fig. As the data show some consider¬ 
able scatter, we also plot two envelope curves (shown as 
dashed lines in Fig. 13, which define a maximally possi¬ 
ble range for the transition pressure from Pmin to Pmax- 
These two parameters are also given in Table and were 
derived such that they define the minimal range of tran¬ 
sition pressures for which no data points fall outside the 
envelope curves. Due to the above-mentioned hysteresis 
effect and vibrat ions d uring handling in the laboratory 
experiments (see 2.3.2), which result in a too-low pack¬ 
ing density for some of the data points and a compaction 
above RCP for others, respectively, the above conditions 
is strictly obeyed for the microgravity data, but we allow 
the laboratory data to fall outside the envelope curves. 

In Fig. all derived transition pressures Pm and their 
respective uncertainty ranges, Pmin — Pmax (^1 data iden¬ 
tical to those of the first three rows in TableQ, are shown 
as asterisks with error bars as a function ofthe regolith- 
particle radius. The data po int from the omnidire ctional 
compression ei^eriments of Guttler et al. (2009|) is also 
shown in Fig. ® at = 0.75_um (asterisk without error 
bar). The solia line in Fig. [^corresponds to Eq. [^ i.e. 

Pm oc . The dotted line is the best fit to the data as¬ 
suming a power law and minimizing the chi-square error 
statistics of the logarithmic data points. 

We think that this deviation in slope can be attributed 
to the rather large errors in our measurements, but is not 
proof that the model behavior (Eq. ^ is wrong. 

The rough agreement between model (Eq. ™ and data 
indicates that rolling is the major effect at tte particle 
level during compression and shows that it can describe 
our measurements reasonably well. Although the Pm was 
not used to fit the mo del to our data , the tr ansition pres¬ 
sure Pm measured by Guttler et al. (2009) for /rm-sized 
monomer grains is alfo reasonably well represented by 
the model. 


3.2. Regolith consisting of dust agglomerates 

At low pressures, the total filling factor of the regolith 
consisting of dust agglomerates should be 


^total ^aggl^RSS: 


( 12 ) 


with $aggi=0.35 (Weidling et al. 2012) and ^rss being 
the internal filling tactor ot tlie iUU-250 /rm-sized dust 
agglomerates and of the regolith super structure, respec¬ 
tively. The force on a single dust agglomerate inside the 
regolith can be calculated from the pressure E within the 
regolith, analogous to Eq. iby 




E TT 

agg 

^RSS 


(13) 


We assume that the friction against rolling-force for an 
agglomerate compared to solid particles of the same size 
is smaller by a factor of $aggi and the regolith super 
structure was formed by random ballistic deposition, 
which yields a filling factor of $j^5s=0.15. We can then 
calculate the minimum pressure at which rolling occurs 


Proll — 


Froll ^RBD 

^ ^aggl ^aggl ■ 


(14) 


We get Proll « 10 ^ Pa and this is the pressure at 
which restructuring of regolith consisting of agglomer- 


10 ® 


10 '' 


" 10 ^ 
CL 

E 

Q. 

10 ° 


10 “^ 


lO"'' 10"® 10"® lO"" 10"® 

Particle Radius [m] 

Fig. 6.— The transition pressure pm from Eq. [^as a function of 
the radius of the spherical regolith particles. The solid line shows 
the model curve, given by Eq. The data point for the smallest 
particle size stems from the compression measurements o f|Gu^ler| 
|et al.| l |2009[ | . All other data points are derived by fitting 
our measurements. The error bars denote the envelope curves in 
Fig. El and are given by pmiFi and pm<^x. 



ates should start. This is consistent with our measure¬ 
ments shown in Fig. 0 (right panel) and the analytical 
approximation shown in Fig. 0 (solid lines). 

3.3. Limits of our Model 

Our model is based on measurements on regolith for 
spherical grains with radii between = 7 x 10“^ m and 
Cp = 5 X 10“^ m and pressures in between E = 10“^ Pa 
and E = 10® Pa. We apply our model to celestial objects 
for which the representative regolith-grain sizes modeled 


to be between ^ = 4.4 x 10 


-5 


m and = 4.2 x 10 


-2 


m (Gundlach et al. 2011). The pressures at the sec¬ 
ond regolitli^^tiele layer in these objects are between 
E = 10“^ Pa and E = 10 Pa and can reach values of 
E = 10® Pa deep inside the regolith. Thus, the parti¬ 
cle sizes have to be extrapolated by only one order of 
magnitude using the results of our experiments and the 
pressures are well within our measurement range. The 
fact that we only measured at a minimum i^-level of 0.06 
ms“^ and, e.g., the asteroid 1996FG3 possesses a g-level 
of 3.3x10“'*^ ms“^ at its surface is not a severe restric¬ 
tion of applicability, because only the pressure inside the 
regolith is relevant. 

However, regolith particles on planetary objects are 
not mono-disperse spherical grains but are irregular in 
shape and possess some size distribution, a s can be seen 
in the Apollo samples (Heiken et al. 1991). To test the 
influence of grain irregularity and poly-dispersity on the 
mechanical properties of the regolith, we previously al¬ 
ready did unidirectional compression measurement s with 


monodispe rse and polydisperse irregular grains ( Blum 
.| 2006|) with a size range of 0.1p,m - 10/rm (see Kothe 
l.l|2(J13l their Fi 


et al 

et al 


Figure 3). We found that the filling 
tactor tor liBD drops from $i = 0.15 for monodisperse 
spherical grains to $1 = 0.11 and $1 = 0.07 if mono¬ 
disperse irregular and poly-disperse irregular grains are 
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TABLE 1 

Parameters of the fit function (Eq. for the four particle types used in this study (rows 1-4), extrapolations to 

MM-SIZED and CM-SIZED AGGLOMERAT E LAYERS (ROWS 5-6) AS WELL AS FOR OTHER COMPRESSION STUDIES (ROWS 7-8). PLEASE NOTE THAT 
(/)1 FROM THe|MACHH ET AL.|(|2013|) DATA WAS INCREASED TO FIT TO OUR MEASUREMENTS. ROWS 9-11 SHOW THE PARAMETERS FOR THE 
MODEL CURVES WITH CONSTANT A (SOLID CURVE IN Ficj^, WITH A OC r®'® (DOTTED CURVE IN FlC. AND WITH A OC r”®'® (DASHED 

CURVE IN Fig. All pressures are given in Pa. 


Particle type 


<t>l <t>2 


Pm 


10 spheres 
60 /rm spheres 
1 mm spheres 
100-250 fim agglomerates 
1 mm agglomerates 


0.15 

0.15 

0.15 

0.05 

0.05 

0.05 


0.64 

0.64 

0.64 

0.37 

0.55 

0.55 


5x10® 

0.35 

0.3 

0.6 

6.1 X 10' 
4.7 X 10' 


Pmi: 


Pmax 


Power 


40 

0.02 

0.005 

0.06 


1.3x10® 

16 

10 

6 


0.58 

0.58 

0.58 

0.58 

0.58 

0.58 


Guttler et al. (|2009ll 

Machii et al. l2U13l) 

0.12 

0.58 

13 X 10® 

0.58 

0.35 

0.55 

4 X 10® 

0.58 

Solid curve in Fig. [5j 

0.15 

0.64 


0.58 

Dotted curve in Fig. 

1 0.15 

0.64 


0.58 

Dashed curve in Fig. 

[ 5 ] 0.15 

0.64 


0.58 


357 Pa 
357 Pa 
357 Pa 


0 

-1-0.5 

-0.5 


used (Blum et al. 2006). The maximum filling factor 
at high unidirectional compressions are $2 = 0.33 for 
mono-disperse spherical and irregular grains and drop 


to $2 = 0- 20 for poly-disperse irregular particles (Blum 
et al.|2006 ). Mind that these values are much lower than 
the KCB value, because under unidirectional compres¬ 
sion the samples creep sideways. Although not explic¬ 
itly given in our previous paper, the values for pm are 
almost identical for the three types of grains and read 
Pm = 5 X 10^ Pa, pm = 2x 10^ Pa, and pm = 4 x 10^ Pa 
for the monodisperse spherical, the mono-disperse irreg¬ 
ular and the poly-disper se irregula r grain s, respectively. 

The measurements of Yu et al. (19971 show that ir¬ 
regular particles behave like s phCT ical grains that are a 
factor of 5 smaller (see Section |3.1l). Applying Eq. we 
find a corresponding change of Pm by a factor of 8. How¬ 


ever, as a size distribution (as e.g. found by Miyamoto 
et al.|2007 ) reduces this effect on pm by a factor of 2 (see 
above), we therefore conclude that the morphology and 
poly-dispersity of the grains play only a minor role for 
the compression behavior of regolith. 

The stratification measurements of the densities of lu¬ 
nar re golith show externa l packing densities larger than 
RCP (jlleiken et al. 1991). This indicates that polydis- 
perse grains can be compacted beyond this limit and 
show that our model underestimates the packing den¬ 
sities by approximately 12% (discussed later). 

Regolith from particles larger than 100 pm is sensitive 
to vibrations due to the large inertia of the particles com¬ 
pared to their inter-particle van der Waals forces. Thus, 
vibrations caused by impacts on the celestial o bjects can 


incre ase their regolith packing density by 0.1 (|Yu et al. 
1997 their Figure 5). 


We do not think that sintering is an important issue 
for regolith particles, because these grains were produced 
by impact destruction and size selective re-accretion. 
Whether the re-accreted particles are directly formed 
by impact fragmentation or by fragmentation and re¬ 
sintering, is not of relevance. The latter is more realistic 


as the inve stigations in the Lunar Sourcebook (Heiken 


et al. 
terna 


1991) show that the regolith particles have a m- 
packing density of about 79% from void internal 


spaces. During the formation of comets in the early so¬ 
lar nebula, the temperatures were presumably too small 
for relevant sintering processes, so that sintering will not 
have any influence for the stratification of cometary sur- 


face regol i th. Ho wever, laboratory measurements of |Pat-| 
El et al. (2009) show that sintering during an orbital 
passage of the comet can increase the stability of the 
regolith in the uppermost few centimeters. This means 
that our model may not be applicable to the uppermost 
centimeters below the comet surface. However, as dur¬ 
ing each orbit a comet can lose surface material of a few 
meter depth, it is not clear whether heat-flow induced 
sintering is of relevance. 

4. APPLICATIONS TO PLANETARY SCIENCES 
4.1. Rocky objects 

In this Section, we will derive the filling factor of the 
regolith on atmosphere-free planetary bodies as a func¬ 
tion of depth below the surface. With this result, we 
will calculate the mechanical yield of the regolith under 
a mass of 100 kg resting on the surface on feet with a 
total area of 0.07 m^. We will consider the effect of am¬ 
bient gravity only and will disregard the compaction of 
the regolith by high-velocity impacts. In that course, we 
will compare our results with measurements taken by the 
Apollo 15-17 core tubes and drill cores on the Moon. 

We again start with the calculation of the pressure 
inside the regolith of an arbitrary planetary body as a 
function of depth h below the surface, T,{h). Because 
the function S(/i) is non algebraic, we obtain its inverse 
h(T,) by integrating the ansatz 


dh = 


dY.{h) 


ptg^E) ’ 
using $(E) from Eq. and get 
E 

$ 2 ESp^ 


(15) 


h = 


+ (*--*.) X (16) 

}■ 


H2F1 




A, 1,1 +A, 


$1 - 2^2 


with H2F1 being the Gaussian hypergeometric function. 
Inserting the inverse function of Eq. 


S($) =pm 


$2 - $1 


A In 10 


(17) 
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Fig. 7 . — The filling factor of the regolith of various Solar-System 
bodies as a function of hydrostatic pressure. The curves from left to 
right Phobos (2.2 mm), asteroid 1996FG3, asteroid Steins, asteroid 
Dodona, asteroid Vesta, the Moon and planet Mercury respectively. 

we get 


h = 


Pm 


$1 - $ 


gpbi^i^2 - 2$^) V® - 

— 2$2 (^2 — ^ l ) ^ 


A 


(18) 


H2F1 


A, 1, 1 + A, $2Pm 




(j) _ $ 


-Pm 


]■ 


The filling factor as a function of pressure and depth be¬ 
low the surface of the regolith is shown in Figs. [7] and 
respectively, for a variety of planetary bodies with dif¬ 
ferent gravitational accelera tions q. The particle sizes o f 
the regolith were taken from Gundlach and Blum| (|2013|) . 
For the graphs shown in both Figures, we used =1X15 
and ^2 = 0.64. is calculated from the particle diame¬ 
ters of the regolith, using Eqs. [^nd [^ (see Fig. |^. The 
graph in Fig. |^is the inverse of Fq. 

Because all reasonable loads, like, e.g. a 100 kg ob¬ 
ject on the surface with a foot-area of 0.07 m^, lead to 
pressures that are orders of magnitude larger than the 
pressures needed to compact the regolith to RCP, we did 
the ansatz for the total yield of the regolith 






^yield — 


di7- 


^RCP 


(19) 


Because Eq. is monotonically increasing and contin- 
uou s, w e achieve the primitive of the inverse function of 
Eq. 19 (second integral) with the substitution H = /i($) 
and partial integration, i.e. 


fffo. 


<^{H)dH = i7oo^(i7oo) - 0$(0) (20) 




Executing th e int egration of the first term in Eq. 19 and 
inserting Eq. 20 results in 


/lyield = / h($)d$, 

J <J»o 


( 21 ) 


which we solved numerically. The above term converges 
very slowly with $oo —t ^rcp = 0.64. This implies that 



Fig. 8 .— The filling factor as a function of depth of the regolith 
(with particle diameters given in parentheses) of (from top to bot¬ 
tom) Mercury (44 /im), the Moon (96 /i.m), Dodona (0.6 mm), 
Vesta (108 /J-m), Phobos (2.2 mm), Steins (1.3 mm), 1996FG3 
(2.0 mm), respectively. 


^yieid depends on the regolith layer depth above a solid 
planetary interior. However, granular regolith has the 
characteristic to distribute point loads on its surface in a 
cone-shaped manner to a larger and larger area beneath. 
These cones hav e typically a frictio n angle between 35 
and 50 degrees (Wu and Sun 2008). This means that 
after a given depth inside a granular media, the force is 
being distributed over a large area. If the correspond¬ 
ing pressure reaches the elastic limit of the regolith, the 
compaction terminates. Here, we assume that the fric¬ 
tion angle is 35 degrees and that the elastic limit of the 
regolith is given by Eq. In Fig. we show the yield 
as a function of the layer aepth (using Eq. 211 for differ¬ 
ent planetary bodies. We plot the curves only for depths 
for which the pressure is above the elastic limit of the 
regolith. As the depth of the regolith is basically un¬ 
known, the maximum possible yield is given by the right 
end of each respective curve. However, if the regolith 
layer is shallower than that and is followed by bare rock 
beneath, the respective yields are smaller and can be de¬ 
rived from the curves in Fig. One can recognize from 
Fig. [^that even for the deepest possible regolith depths, 
the static yield of the regolith under the weight of the 
assumed spacecraft or astronaut is on the order of a few 
decimeters. 

To compare our results with measurements taken on 
the M oon with Apo l lo 15- 17 core tubes given in Table 
9.5 of H^iken et al. (1991), we average o ur lunar data 
of Fig over the lepth rang es given in IHeiken et liLl 
(|1991 ). The porosity values of IHeiken et al.l (|199ip in¬ 
clude the internal porosity of the lunar regolith grains. 
This in ternal packing dens ity was measured and reads 
^0.79 (Heiken et al. 1991). In Table [2| we present our 
data on the tilling factor of the lunar regolith together 
with the external packing density, corrected for the in- 
ternal packing den sity of the grains (both taken from 
Heiken et aHl99T I. The lunar surface regolith obviously 


is packed sightly denser than predicted by our model. 
This is because of the size distribution and irregularity 
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spacecraft) on different celestial bodies. The graphs give the yields 
(left to right) on 1996FG3, Phobos, Steins (almost identical to 
Phobos), Vesta, Dodona, the Moon and Mercury, respectively, as 
a function of the depth of the regolith layer. Please note that due 
to the friction angle, the yield will not increase further at a regolith 
depth at which the pressure is reduced to the elastic limit. The 
curves terminate at this point. 


TABLE 2 

Comparison of our predicted filling factor (center 
column) of the lunar regolith (averaged over the depth 

RANGE SHOWN IN THE LEFT COLUMN) WITH MEA SURED EXTERNA L 
FILLING FAGTORS F ROM A pOLLO PUBLISHED IN IHeIKEN ET AL.I 
(|1991[l (right column). 


Depth Range (m) 

Our Model 

Measured 

0 - 0.15 

0.57 

0.60 

0 - 0.30 

0.59 

0.64 

0.30 - 0.60 

0.62 

0.70 

0 - 0.60 

0.60 

0.68 


of the lunar regolith par ticles for which packing den sities 
u p to 70% are possib le (jPesmond and Weeks||2013 1. 


Magri et al. (20011 estimated the porosity of the first 
meter ot regolith on asteroids using ground-based radar 
measurements. However, their data possess huge errors 
and it is not possible to find a correlation between fill¬ 
ing factor and size of the asteroids in their data. If we 
average over all their findings, we get a filling factor of 
0.49 ± 0.14. If we assume that the regolith grains of 
the asteroids possess the same internal packing density 
as measured for the lunar regolith particles, the average 
packing density is 0.62 ± 0.14, which matches with our 
model that predicts that at a few decimeters depth the 
regolith is at RCP. This estimations are, however, not a 
strong confirmation of our model, because the averaging 
over objects of different sizes and taxonomic types might 
not be appropriate. 


4.2. Comet 67P/Churyumov-Gerasimenko 


According to the model of Skorov and Blum (2012), 
comet nuclei consist of dust and ice agglomerates, each 
consisting of /rm-sized solid dust and ice grains. To cal¬ 


culate for this model the pressure, depth and yield-load 
dependencies of the fillin g fac tor on a comet, we follow 
the calculation in Section 4.1 Here, we assume that ice 
agglomerates have a transition pressure that is a factor 


of 10 larger than that of our silicate agglomerates (|Gund- 
lach et aL]|2011 ). We calculated the comet-regolit' 


n com- 

pression tor agglomerates of 1 cm, 1 mm, 180 ^m in size, 
respectively, the latter being the dust-aggregate size in 
our experiments. 

Particles that are in contact possess a contact surface 
whose size is proportional to their radius squared. Dur¬ 
ing rolling, additional contact surface is added in the di¬ 
rection of rolling and removed in the opposite direction. 
However, this process of creation and removal of contact 
surface is asymmetric, because the contact forces pull the 
contacting surfaces outward before they break. This dif¬ 
ference causes a to rque on the particle, which is known as 


rollin g friction (see Krijt et al. 2014 Dominik and Tielens 


19951. The radius dependence ot this mechanism is found 


by Krijt et al. (2014) to be given by ra. Agglomerates 
in contact possess a contact area for which the number 
of monomer particles is proportional to the agglomer¬ 
ate radius squared. During rolling, new particle-particle 
contacts are formed in the direction of rolling and broken 
in the opposite direction, because the monomer particles 
are also pulled outward. Because of this analogy be¬ 
tween solid-solid and agglomerate-agglomerate contacts, 
we assume that the agglomerate-radius dependence of 
this mechanism follows the same power law as for com¬ 
pact particles. The parameters used in our calculations 
are also given in Table 

In the following, we apply our calculations to comet 
67P/Churyumov-Gerasimenko, which is being visited by 
the Rosetta spacecraft in 2014/2015. As the Rosetta- 
mission data concerning the g-levels have not been pub¬ 
lished at the submission date of this paper, we calcu¬ 
lated them using the published mass, average mass den- 
sity and volume of comet 67P/Churyumov-Gerasimenko 
(Sierks et al. 2015) of m = 10^^ kg, p = 470 kg m“^ 
and P = 25 km'^. Thus, we get for the average radius of 
i? = 1.8 km a surface acceleration of g = 2 x 10“'* m s“^. 

We calculated the properties of cometary material for 
two extreme cases, namely that he comet consists of dust 
agglomerates or ice agglomerates only, and one average 
case that the comet consists of a 50%-50% (in volume) 
mixture of dust and ice agglomerates, as proposed by 


Skorov and Blum (2012). We further assume that the 
dust and ice agglomerates are arranged in random or 
alternate order so that the compressive behavior can be 
averaged in the following way 


4^ice-t-dust (p)— 2 (*^ice(p) T 4^dust(p)) 

Yieldice-Fdust (p) = ^ (Yieldice(p) + Yielddust(p)) ■ (23) 

Fig. shows the filling factor as a function of pressure 
for comet-nucleus material, assuming that the comet con¬ 
sists of agglomerates of 1 cm size (first group of 3 curves), 
1 mm size (second group of 3 curves) and 180 pm size 
(third group of 3 curves). Within these groups, the ag¬ 
glomerates consisting of dust agglomerates are always 
represented by the left curve, the ice-and-dust mixture 
by the center curve, and the ice agglomerates by the right 
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curve, respectively. 

Regolith consisting of agglomerates has a smaller total 
RBD filling factor compared to regolith from solid grains, 
due to the intrinsic porosity of the agglomerates (see 
Section 13.21). Therefore, they can be compacted further 
than RGB of the regolith super structure by destruction 
or deformation of the spherical agglomerates and subse¬ 
quent compaction to a global RGB structure. Its strength 
against compression is higher than that of a regolith con¬ 


sistin g of compact grains of the same size (Machii et al. 


2013). The compression curve thus consists of two parts. 


ihe first compaction stage is due to the reorientation of 
the agglomerates by rolling until RGB of the s uper struc¬ 
ture is reached (our measurements in Section 2.3). The 
data for the second compact ion stage are based on the 
compression measurement of Machii et al. (2013), which 
follow the compression of the agglomerates themselves. 
This compression takes places at pressures larger than 
10^ Ba, because the /rm-sized grains (i.e. the constituent 
grains of the agglomerates) d etermine the compr ession 
strength. The original data of Machii et al. (2013) start 
at a filling factor of 0.24, whereas the maximum com¬ 
pression of our measurements ends at a filling factor of 
0.37. This is due to the different preparation methods of 
the agglomerates. The more spherical and monodisperse 
the a gglomerates are, th e smaller is their filling f actor at 
RGB (jPonev et al. |2004[ ) . We used the results by lMachli 


et al. (2013) by increasing their initial filling factor to our 


end value of 6.37 (see Table and assumimg that the 
pressure required for the compression of icy agglomerates 
to a given filling factor is 10 times higher than for dusty 
ones. For compressions where the agglomerates are in 
RGB, the total filling factor as a function of pressure is 
independent of the agglomerat es size. This means that 
the curve corresponding to the Machii et al. (2013) data 
(slopes at high pressures) is valid tor all three agglom¬ 
erate sizes. An explanation for this behavior is given in 
the Appendix. 

Mind that the pressures at the second compression 
stage are much higher than the internal pressures on 
small Solar-System objects. Such pressures can occur 
on a comet only during high-velocity impacts. For im- 
p act pressures in the GBa range we recommend the work 


of Beitz et al. 
4’he high ro 


(2013). 


ling and sticking forces of water ice lead to 
a comet nucleus of low filling factor; whereas the regions 
close to the comet surface presumably consist of dust 
agglomerates only, due to the sublimation of water ice, 
and will thus possess a higher filling factor. A mixture 
of dust and ice agglomerates possibly results in a filling 
factor in between the extreme values (see Fig. 10). A 
comet nucleus is assumed to possess a few centimeter 
thick layer of pure dust agglomerates, which is followed 
by a mixture of dust and ice agglomerates. Thus, the 
filling factor will follow the respective left curve of one 
of the three groups in Fig. (10) in the upper centimeters 
of the surface, followed by the center curve for a 50%- 
50% (in volume) mixture of dust and ice. For a higher 
ice content, the filling factor will be closer to the pure 
ice curve, which is the right-most curve in the respective 
triples. For a lower ice abundance, the filling factor will 
be c lose r to the left curve. The three triple-curves in 
Fig. IT represent dust and ice agglomerate sizes of 1 cm, 
1 mm, and 180 /rm, respectively. 



Fig. 10.— Filling factor as a function of pressure of cometary 
material. The first group of three curves represents regolith from 
cm-sized agglomerates, the second group of three curves represents 
regolith from mm-sized agglomerates, and the third group of three 
curves represents regolith from 180 pm-sized agglomerates, respec- 
tivel y. The high-pressure behavior of the material is based on data 
from|Machii et al.|||2013|l. These curves are independent of the ag- 
glomerate size ot the regolith. Within these groups, the left curve 
is regolith from dusty material, the middle curve is a 50%-50% (in 
volume) mixture of agglomerates from dust and water ice, and the 
right curve is regolith from water ice. Please note that the curve of 
the ice particles of the mm-sized agglomerates is almost identical 
with the curve of the dust particles of the 180 pm sized particles. 


Applied to comet 67P/Churyumov-Gersimenko, we get 
a filling factor as a function of depth as shown in Fig. 
[g The respective triple curves are, those for agglomer¬ 
ate sizes of 1 cm (dotted lines), 1 mm (dashed lines), and 
180 /rm solid lines. The respective upper curve represents 
agglomerates consisting of dust, the center curve repre¬ 
sents 50%-50% (in volume) mixtures of dust and ice, and 
the lower curve stands for pure ice agglomerates, respec¬ 
tively. All curves start from a layer depth of one particle 
diameter. 

To apply our results to the anticipated landing of the 
Rosetta lander on comet 67F/Ghuryumov-Gerasimenko, 
we calculated how much the surface of the comet nucleus 
will yield at touchdown. For the lander, we assume a 
mass of 100 kg and a total landing-feet area of 0.07 m^. 
In case the spacecraft rests on the surface it will only 
render the agglomerate super structure to RGB and will 
not considerably sink in under its own weight, because 
the pressure below its feet is only ~ 0.4 Ba and, thus, 
much lower than the pressure needed to compact the re¬ 
golith (see Fig. 10). However, the impact pressure will 
be higher than the static weight of the sp acecraft . _ 

Instead of using an impact model (e.g. Melosh 1989), 
we estimate the impact pressure by assuming that the 
impact velocity will be decelerated steadily by the land¬ 
ing gear on a length of 10 cm, which is foreseen for 
the Rosetta lander. This dynamic pressure is treated 
as a st atic pressure on the surface. According to |Melosh| 
(1989), dynamic compression of our regolith would re- 
quire a 20 times stronger deceleration and does therefore 
not appear. An issue of the static ansatz is the small 
deceleration time (0.4 s), because of t he low sound spee d 
in the regolith of about 20 ms“^ (see Beitz et al.||2013). 
This limits the depth to which the regolith can be com¬ 
pressed to 8 m. But due to the friction angle, the elastic 
limit is already reached for 180 ^m-sized aggregates at 
3 m, 7 m and 10 m for pure ice, dust-ice mixture and 
pure dust, respectively, assuming an impact velocity of 
1 m s“^. For the mm and cm-sized agglomerates the elas¬ 
tic limit is reached at a depth larger than 100 m by the 
gravitational acceleration of the comet only. For these 
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Fig. 11.— The filling factor as a function of depth below the sur¬ 
face of Comet 67P/Churyumov-Gerasimenko. The group of three 
dotted curves represents particle layers consisting of 1 cm-sized ag¬ 
glomerates, the group of dashed graphs represents particle layers 
consisting of 1 mm-sized agglomerates, the group of solid graphs 
represents particle layers consisting of 180 ^m-sized agglomerates. 
Within these groups, the upper curve represents dusty material, 
the middle curve a 50%-50% (in volume) mixture of agglomerates 
from dust and water ice, and the lower curve ice agglomerates, 
respectively. All curves start from a depth of one agglomerate di¬ 
ameter. 

agglomerate sizes the yield is independent of the impact 
velocity. _ 

In Fig. the solid group of three lines show the yield 
of a spacecraft landing on regolith consisting of 180 /im- 
sized aggregates as a function of impact velocity. The 
dashed and dash-dotted group of three lines show the 
impact penetration for mm-sized agglomerates and cm 
-sized agglomerates. Within these groups of three lines 
the uppermost line corresponds to ice agglomerates the 
center line corresponds to a 50%-50% (in volume) mix¬ 
ture of ice and dust agglomerates and the lowermost line 
corresponds to dust agglomerates.. 

At low velocities (fiimp <20 cm s“^) and a surface of 
180 ^m-sized aggregates , the yield is dominated by the 
gravity of the comet. For larger velocities the impact 
velocity gets important and dominates the yield at ve¬ 
locities larger than 1 m s“^. In case the surface consists 
of agglomerates smaller than a few millimeters in diam¬ 
eter the yield is always dominated by the gravity of the 
comet. 

The penetration depth is higher for the ice agglomer¬ 
ates, because their higher rolling force prevented an RCP 
packing down to a depth of a few decimeters so that the 
icy material can potentially be compacted more than a 
pure dust-aggregate surface. However, the total penetra¬ 
tion is less than 20 cm even for impact speeds as high as 
1 m s“^. 



Fig. 12.— The yield of the Rosetta lander (assuming a mass of 
100 kg and a total landing-feet area 0.07m^) whose landing gear 
decelerates the impact within a length of 10 cm on an agglomerate- 
layer material consisting of 180 ^m-sized agglomerates (solid lines), 
1 mm-sized agglomerates (dashed lines) 1 cm sized agglomerates 
(dashed-dotted lines), as a function of impact velocity. Within 
a group of this lines the uppermost curve is in case the surface 
consists of ice agglomerates, the lowermost curve in case of dust 
agglomerates and the center curve in case of a 50%-50% (in volume) 
mixture of dust and ice agglomerates. 


the filling factor depends on the ambient hydrostatic 
pressure as well as on grain size and morphology. Using 
an analytical description for the compression curve, we 
developed an analytical model to predict the stratifi¬ 
cation of regolith-covered dusty and icy Solar-System 
bodies. We compared our findings of the stratification 
of the surface regolith of the Moon with Apollo 15-17 
core tubes and drill cores and found only a slightly 
denser packing than predicted in our model. In a second 
step, we calculated the mechanical yield of a spacecraft 
resting on the surface of such a body. We found that 
a spacecraft with reasonable mass will never penetrate 
into the regolith of any Solar-System body by more than 
60 cm and a comet by more than 25 cm. 
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5. SUMMARY 

From measurements of the static compaction of 
regolith analog at different g-levels on-board a parabolic 
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APPENDIX 


The calc ulation in this Sectio n shows that the high-pressure behavior of the filling factor of dust agglomerates, as 
studie d by Machii e t al. (2013) and represented in Fig. 10, is independent of the size of the agglomerates. We as 


well as Machii et al. 12013) used in the experiments agglomerates with a filling factor of 0.35. These agglomerates are 


organiz ed m RCF in a regolith with a given minimal pressure (which is depending on the agglomerate size, as shown 
in Fig. 10). Because the agglomerat es can fill the void spaces in between each other under further compression, the 
compression curve follows the one by Blum and Schrapler (2004) for unidirectional compression. The pressure above 


which the grains inside the agglomerates flow is p > pf = 10°Pa Blum and Schrapler (see 2004). 

Therefore, the area of the contact circle between two agglomerates, Ac, depends on the force applied to the contact 
sphere. 


Pf 


(Al) 


The dependency of the contact area on the indentation h of the deformed agglomerates with radius r is approximately 
given by the cut-face of a sphere 

Ac = {2hr-h^). (A2) 


Applying an indentation that is proportional to the radius of the agglomerates, i.e. h/r = const, we get Ac oc r^. The 
force on the agglomerates inside the body depends on the local pressure p through 


pr^TT 

^RBD 


(A3) 


(see also Eq. [^. With the above condition, we see that the pressure acting at the contact area of an agglomerate is 
independent m the radius of the agglomerates, because both the force on the agglomerate at a given pressure inside 
the agglomerate and the contact area at a given force are proportional to r^. This also means that a given pressure 
leads to a compaction of the regolith independent of the agglomerate size. This conclusion is, however, only valid for 
pressures high enough to compact the agglomerates, i.e. for p > py. 
































